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I A new parameterization for an effective non-linear Lagrangian den- 

sity of relativistic mean field (RMF) theory is proposed, which is able 
, to provide an excellent description not only for the properties of sta- 

ble nuclei but also for those far from the valley of beta-stability. In 
addition recently measured superdeformed mimima in the Hg-region 
^ . are reproduced with high accuracy. 

P5 ! Relativistic Mean Field (RMF) |l|] theory has recently gained considerable 

success in describing various facets of nuclear structure properties. With a 
very limited number of parameters, RMF theory is able to give a quantitative 
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Abstract 



;_j ■ description of ground state properties of spherical and deformed nuclei 



at and away from the stability line. Recently it has been shown that RMF 
theory is successful in reproducing the anomalous kink in the isotope shifts 
of Pb nuclei and a first-ever microscopic description of anomalous isotopic 
shifts in Sr and Kr chains P] has been provided. Such an anomalous behavior 
is a generic feature of deformed nuclei which include almost all isotopic chains 
in the rare-earth region P] where RMF theory has been shown to have a 
remarkable success. Moreover good agreement with experimental data has 
been found recently for collective excitations such as giant resonances 
and for twin bands in rotating superdeformed nuclei 0. It is also noted that 
cranked RMF theory provides an excellent description of of superdeformed 
rotational bands in the A=140-150 region p[, in the Sr region |10l and in the 



Hg region |11 



The starting point of RMF theory is a standard Lagrangian density 
C = {-f{id - g^u - Qppf -eA) - m- g„o) 
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-\r,uR^'' + \my - (1) 

which contains nucleons if) with mass m, cr-, a;-, p-mesons, the electromag- 
netic field and non-linear self-interactions of the a-field, 

U{a) = Knla^ + ^g^a' + \g^a\ (2) 

The Lagrangian parameters are usually obtained by a fitting procedure 
to some bulk properties of a set of spherical nuclei |jl2|. Among the exist- 
ing parameterizations the most frequently used are NLl |jl3[, NL-SH [|n 



and the parameter set PL-40 |jT5|, which has been proved to provide reason- 
able fission barriers. NLl and NL-SH sets give good results in most of the 
cases. Along the beta stability line NLl gives excellent results for binding 
energies and charge radii, in addition it provides an excellent description of 
the superdeformed bands p, |ll|. However, in going away from the stability 
line the results are less satisfactory. This can be partly attributed to the 
large asymmetry energy J ~44 MeV predicted by this force. In addition 
the calculated neutron skin thickness shows systematic deviations from the 



experimental values for the set NLl [16|. In the parameter set NL-SH this 



problem was treated in a better way and improved isovector properties have 
been obtained with an asymmetry energy of J ~ 36 MeV. Moreover NL- 
SH seems to describe the deformation properties in a better way than NLl. 
However, the NL-SH parameterization produces a slight over-binding along 
the line of beta-stability and in addition it fails to reproduce successfully the 
superdeformed minima in Hg-isotopes in constraint calculations for the en- 
ergy landscape. A remarkable difference between the two parameterizations 
are the quite different values predicted for the nuclear matter incompress- 
ibility. NLl predicts a small value (K=212 MeV) while with NL-SH a very 
large value (K=355 MeV) is obtained. Both forces fail to reproduce the 
experimental values for the isoscalar giant monopole resonances for Pb and 
Zr nuclei. The NLl parameterization underestimates the empirical data by 
about 2 MeV while NL-SH overestimates it by about 2 MeV. 

The aim of the present investigation is to provide a new improved set 
of Lagrangian parameters, which to some extend cures the deficiencies of 
the existing parameterizations. For this reason a multi-parameter fit was 



performed in the the same way as with the other parameterizations [|12|, [14 
The nucleon mass was fixed to 939 MeV. The Lagrangian parameters are 
the meson masses m^j, m^, mp, the corresponding coupling constants go-, g^, 
gp and the parameters g2, gs of the non-linear potential U(cr). Apart from 
the mass of the p meson which was fixed to the empirical value (763 MeV) 
all the others were taken as free parameters. The nuclear properties fitted 
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are the charge radii, the binding energies, and the available neutron radii 
of several spherical nuclei. The experimental input for finite nuclei used in 
the fitting procedure is shown in Table 1 in parentheses. We recall that 
for the determination of NL-SH parameters six nuclei were used in the fit, 
namely ^^O, ^^Ca, ^^Zr, ^^^Sn, ^"^^Sn and ^^^Ph while for NLl ^^Ca and 
^^Ni were also taken into account. It is noted that for NLl the experimental 
information used was the total binding energies, the diffraction radii and 
the surface thickness. For NL-SH charge radii and neutron radii were used 
instead of the diffraction radii and the surface thickness. In the present work 
the number of nuclei used in the fit was increased to ten. In order to take 
into account a larger variation in isospin, in addition to the eight nuclei used 
for NLl the doubly closed shell nucleus ^^^Sn as well as the heavier lead 
isotope ^^^Pb were also included in the fit. The experimental values for the 
total binding energies were taken from the experimental mass tables [IT^, the 



charge radii from Ref. |T^. The available neutron radii are from Ref. [|1£ 
In the case of open shell nuclei pairing was considered in the BCS formalism. 
The gap parameters A„(p) were determined from the observed odd-even mass 
differences [0. Specifically, for ^^Ni, A„=1.4 MeV, for ^^Zr Ap= 1.12 MeV, 
for the two Sn isotopes (A=116,124) the values are 1.17 and 1.32 MeV 
respectively and finally for '^^'^Ph An=0.7 MeV. The binding energies and 
charge radii were taken within an accuracy of 0.1% and 0.2% respectively. 
For the neutron radii, however, due to existing uncertainties the experimental 
error taken into account was 2%. In addition in the fitting procedure some 
nuclear matter properties were also considered. As "experimental input" 
the following values were used: E/A= -16.0 MeV (5%), p = 0.153 (fm~^ 
(10%), K = 250 (MeV) (10%) J = 33 MeV (10%). The values in parentheses 
correspond to the error-bars used in the fit. 

In Table 1 we list the predictions of NL3 for the ground state properties 
of the nuclei used in the fit. It is seen that they are in very good agreement 
with the empirical values. 

In Table 2 we show the values for the new parameter set. Adopting 



the convention introduced by P.-G. Reinhard [|T2], |T3|, [T^ for the non-linear 
parameterizations the set is named NL3. 

In the same table we give nuclear matter properties calculated with NL3. 
The saturation density p has the value 0.1483 fm~^. The effective mass m*/m 
was found 0.6. It is the same as for NL-SH and slightly higher than for NLl. 
The nuclear matter incompressibility has the value K =271.8 MeV. It is 
therefore somewhere in the middle between the values predicted by NLl and 
NL-SH. Finally the asymmetry energy J is 37.4 MeV. It is closer to that of 
NL-SH and much smaller than that of NLl. 

In the following we present some applications of the new parameter set 
NL3 using the various RMF codes of the Munich group. We performed de- 
tailed calculations for the chain of Sn isotopes with the spherical Relativistic 
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Hartree Bogoliubov (RHB) code discussed in Ref. [^. In Fig. 1 we show 
the isotopic dependence of the deviation of the theoretical mass calculated in 
RMF theory from the experimental values for Sn nuclei. The theoretical 
results were obtained using the parameter sets NLl, NL-SH and NL3. It is 
seen that all parameterizations give a very good description of the experi- 
mental masses. It is also seen, however, that the new force NL3 is able to 
provide improved results over the NLl and NL-SH, reducing rms deviation 
of the masses. 

Axially symmetric calculations have been performed for some well de- 
formed rare-earth and actinide nuclei using the new Lagrangian parameter- 
ization NL3. Here the pairing correlations are taken into account using the 
BCS formalism. The pairing parameters A„(p) were taken from tables XI 
and XIII of Ref. 0. In Table 2 we give the results of our calculations 
together with the experimental information whenever available. It is seen 
that NL3 gives excellent results for the ground state properties of rare-earth 



and actinide nuclei. The experimental masses [|17| are reproduced within 
an accuracy of a few hundreths of keV. The charge radii are in very good 
agreement with the experiment [|1^]. The deformation properties are also 
in excellent agreement with the empirical values. The absolute values of 



the empirical P2 were obtained from the compilation of Raman et al. [^T 
The experimental data for the hexadecupole moments of rare-earth nuclei 
are from a very recent compilation by Lobner p2|. Finally the experimental 
data for the proton quadrupole moments were taken from tables XII and 
XIV or Ref.. 0. Next we report some preliminary results for the Giant 
monopole breathing energies of ^osp^ g^j^^ ^^Zr nuclei obtained from gener- 
ator coordinate calculations based on constrained RMF wave functions. A 
detailed study including in addition dynamic RMF calculations will appear 



in a forthcoming publication In Table 4 we show results of calculations 
using the new parameter set NL3 and compare it with experimental results 
and calculations obtained from the sets NL-SH and NLl. It is seen that NL3 
is able to reproduce nicely the experimental values while the other two forces 
fail, either underestimating (NLl) or overestimating (NL-SH) the experiment 
by almost 2 MeV. This is an indication that NL3 has a correct value for the 
nuclear incompressibility. 

Recently, the excitation energy between the ground state band and the 
superdeformed band in ^^^Hg was measured for the first time Extrap- 



olating to zero angular momentum the superdeformed minimum was found 
to be 6 MeV above the ground state. Performing RMF calculations with the 
parameter set NL3 and mapping the energy surface by a quadratic constraint 
we found the superdeformed minimum at an excitation energy of 5.997 MeV 
above the ground state. A detailed study will be published elsewhere ^5 . 



In conclusion our calculations with the new RMF parameterization NL3 
give very good results in all cases considered so far. It is in excellent agree- 
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ment with experimental nuclear masses, as well as the deformation properties. 
For the first time a RMF parameterization reproduces the isoscalar monopole 
energies of Pb and Zr nuclei. This gives us confidence that NL3 can be used 
successfully in future investigations together with the other parameteriza- 
tions. 

One of the authors (G.A.L) acknowledges support by the European Union 
under the contract TMR-EU/ERB FMBCICT-950216. This work is also sup- 
ported in part by the Bundesministerium fiir Bildung und Forschung under 
the project 06 TM 743 (6). We thank Prof. K.E.G. Lobner for supplying us 
with his compilation on hexadecupole moments prior to its publication. 
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Table 1: The total binding energies, charge radii and neutron radii used in 
the fit (values in parentheses) together with the NL3 predictions. 



Nucleus B.E (MeV) ich (fm) r„ (fm) 



160 


-128.83 (-127.62) 


2.730 


(2.730) 


2.580 




40Ca 


-342.02 (-342.06) 


3.469 


(3.450) 


3.328 


(3.370) 


48Ca 


-415.15 (-416.00) 


3.470 


(3.451) 


3.603 


(3.625) 


58Ni 


-503.15 (-506.50) 


3.740 


(3.769) 


3.740 


(3.700) 


902J. 


-782.63 (-783.90) 


4.287 


(4.258) 


4.306 


(4.289) 




-987.67 (-988.69) 


4.611 


(4.627) 


4.735 


(4.692) 




-1050.18 (-1049.97) 


4.661 


(4.677) 


4.900 


(4.851) 


i32Sn 


-1105.44 (-1102.90) 


4.709 




4.985 




208pb 


-1639.54 (-1636.47) 


5.520 


(5.503) 


5.741 


(5.593) 


214pb 


-1661.62 (-1663.30) 


5.581 


(5.558) 


5.855 





Table 2: Parameters of the effective interaction NL3 in the RMF theory 
together with the nuclear matter properties obtained with this effective force. 



M 


= 939 (MeV) 






rria 


= 508.194 (McV) 


9a 


= 10.217 




= 782.501 (McV) 


9lo 


= 12.868 


nip 


= 763.000 (MeV) 


9p 


= 4.474 


92 - 


= -10.431 (fm-^) 


93 


= -28.885 



Nuclear matter properties 



po 0.1483 fm-=^ 

(E/A)^ 16.299 MeV 

K 271.76 MeV 

J 37 A MeV 

m*/m 0.60 
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Table 3: Total binding energies B.E (in MeV) , charge radii rc (in fm), 
quadrupole deformation parameters (32, proton quadrupole moments Qp (in 
barns) and proton hexadecupole (hp) moments (in barns^) for some deformed 
rare-earth and actinide nuclei with the parameterization NL3. The values in 
parentheses correspond to the empirical data. For details see text 



A 


JD .Hi 




P2 


o 


rip 


152 

O III 














-1294.49 


5.177 


0.301 


5.63 


0.48 


158 ry J 


(-1294.05) 


(5.099) 


(0.306) 


(5.78) 


(0.46(2)) 




1 OQR /in 


0. i ( 


U.o4z 


7 1/1 


n /I s 

U.4o 


162 7-) 


(-1295.90) 


(5.172) 


(0.348) 


(7.36) 


(0.39(9)) 




1 'iOA no 


^ 997 
O.ZZ / 


n "^47 

U.041: / 


7 

i .04 


U.40 


166 


(-1324.11) 


(5.210) 


(0.341) 


(7.36) 


(0.27(10)) 




-lOOi.UO 


979 




7 87 


u.oo 




f-1351 57) 


f5 303) 


(0.342) 


(7 70) 


fO 32fl6)) 




-1406.15 


5.336 


0.328 


7.77 


0.04 




(-1406.60) 


(5.410) 


(0.325) 


(7.58) 


(0.22l°i^) 


232^^^ 














-1766.29 


5.825 


0.251 


9.23 


1.06 


236^ 


(-1766.69) 


(5.790) 


(0.261) 


(9.62) 


(1.22) 




-1790.67 


5.873 


0.275 


10.60 


1.16 


238^ 


(-1790.42) 




(0.282) 


(10.80) 


(1.30) 




-1801.39 


5.892 


0.283 


10.93 


1.07 




(-1801.69) 


(5.854) 


(0.286) 


(11.12) 


(1.38) 



Table 4: Isoscalar giant monopole energies calculated with the effective in- 
teractions NL3, NLl, NL-SH along with the empirical values. 



A 


expt. 


NL3 


NLl 


NL-SH 


208p^ 


13.8 ±0.5 


13.0 


11.0 


15.0 


90zr 


16.2 ±0.5 


16.9 


14.1 


19.5 
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